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ERRATA 


Fig*  1  should  be  captioned:  Flexural  waves  from  a  water  shot  (a)  and  air- coupled  flexural 
waves  from  an  air  shot  (b)  in  smooth  ice  1.1  ft  thick. 


Fig.  2  should  be  captioned:  Dimensionless  phase  and  group  velocity  curves  for  flexural 
waves  in  a  floating  ice  sheet. 

In  Fig.  2,  the  upper  curve  represents  U/v  ,  the  lower  curve,  cl v  . 
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In  Fig.  3,  the  cross  plotted  at  /  -  34.  U  - 
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ABSTRACT 


PART  I 

Experimental  studies  of  the  propagation  of  elastic  waves  on  float¬ 
ing  ice  sheets  were  made  on  the  ice  of  Lake  Superior  and  Lake  Cayuga. 
Elastic  waves  were  produced  by  small  explosive  charges  detonated  at 
various  depths  in  the  water,  within  the  ice,  and  in  the  air.  Seismic 
detectors  consisting  of  a  spread  of  geophones,  a  microphone,  and  a 
hydrophone  recorded  the  resultant  wave  motion  at  various  distances. 

Shots  in  the  water  produced  the  normal  sequence  of  dispersive 
flexural  waves.  For  shots  in  the  air  the  dispersive  flexural  waves 
were  absent  and  a  train  of  constant  frequency  waves  was  observed,  be¬ 
ginning  gradually  at  the  approximate  time  t  -  r/ 2v a  and  culminating 
with  the  arrival  of  the  air  wave  at  the  time  t  =  r/v  (r  is  the 
range  and  va  is  the  speed  of  sound  in  air).  These  waves  were  inter¬ 
preted  as  air-coupled  flexural  waves.  The  frequency  of  the  air- 
induced  waves  is  that  of  flexural  waves  whose  phase  velocity  equals 
the  speed  of  sound  in  air.  The  generated  wave  train  precedes  the  air 
disturbance  since  the  group  velocity  of  flexural  waves  exceeds  the 
phase  velocity,  all  in  accordance  with  classical  theory.  The  fre¬ 
quency  of  the  air-coupled  vibration  is  simply  related  to  the  ice 
thickness. 

The  interpretation  was  supported  by  subsequent  tests  which  con¬ 
sisted  of  shooting  on  shore  and  recording  on  the  ice,  shooting  and  re¬ 
cording  on  ice  sheets  of  varying  thickness,  and  recording  air-coupled 
flexural  vibrations  with  microphones. 

PART  II 

The  exact  theory  for  air-coupled  flexural  waves  in  a  floating  ice 
sheet  is  derived  for  the  case  of  an  impulsive  point  source  situated 
either  in  the  air  or  in  the  water  beneath  the  ice. 
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AIR-COUPLED  FLEXURAL  WAVES  IN  FLOATING  ICE 

PART  I 

EXPERIMENTAL  STUDIES  OF  AIR-COUPLED 
FLEXURAL  WAVES 

INTRODUCTION 

Ekiring  a  period  of  three  weeks  in  February  and  March  1950,  experimental  studies  of  the 
propagation  of  elastic  waves  in  floating  ice  sheets  were  made  on  Lake  Superior  and  Lake 
Cayuga.  These  experiments  were  designed  to  study  comp r ess ion al,  shear,  flexural,  and  lon¬ 
gitudinal  surface  waves  in  the  ice-water  system.  In  the  course  of  the  investigations  an 
air-coupled  flexural  wave  was  discovered.  Sufficient  time  remained  to  firmly  establish  the 
existence  of  this  wave  by  performing  additional  tests. 

The  experimental  data  obtained  for  air-coupled  flexural  waves  and  a  brief  summary  of 
the  theory  are  given  in  Part  I  of  this  report. 

A  detailed  theoretical  investigation  of  these  waves  is  presented  in  Part  II. 


EXPERIMENTAL  PROCEDURE 


The  standard  method  of  seismic  refraction  measurements  was  employed  in  the  experiments 
on  Lake  Superior  and  Lake  Cayuga.  A  spread  consisting  of  8  geophones  and  2  horizontal 
seismometers  was  placed  on  the  ice  surface.  Additional  data  were  occasionally  obtained  from 
a  hydrophone  at  various  depths  in  the  water  and  a  microphone  in  the  air.  Shots  consisted 
of  blasting  caps  or  j  lb  TNT  demolition  blocks  detonated  at  various  depths  in  the  water,  in 
the  ice,  and  at  various  heights  in  the  air.  The  spread  interval  varied  during  the  course 
of  the  experiments,  and  shot  distances  ranged  up  to  7,414  ft.  The  instant  of  the  explosion 
was  transmitted  to  the  recording  truck  by  wire. 

Two  types  of  ice  were  encountered  during  the  tests.  Near  shore  a  rough,  hummocky  ice 
sheet  was  found  in  which  thickness  variations  of  up  to  50%  were  not  uncommon  along  a  pro¬ 
file.  About  a  mile  off  shore  a  smooth  newly  formed  ice  sheet  was  found.  Ice  thickness 
varied  15%  or  less  in  the  smooth  sheet. 


•Manuscript  received  for  publication  18  July  1950. 


10 


FLEXURAL  WAVES 


A  typical  record  of  flexural  waves  for  a  shot  consisting  of  a  blasting  cap  exploded  at 
the  under  surface  of  the  ice  at  a  spread  distance  of  4614-5014  ft  is  shown  in  Fig.  la.  The 
ice  thickness  averaged  1.1  ft  and  the  water  depth  126  ft. 

Flexural  vibrations  in  a  floating  ice  sheet  are  essentially  surface  waves  analogous  to 
gravity  waves  on  water,  Rayleigh  waves  on  the  surface  of  an  elastic  solid,  or  to  flexural 
waves  in  a  rod  or  thin  plate.  Ewing  and  Crary1  made  a  detailed  investigation  of  these 
waves.  They  derived  equations  for  phase  and  group  velocity  as  a  function  of  frequency, 
taking  gravity,  compressibility  of  water,  and  finite  water  depth  into  account.  Their  ob¬ 
servations  of  dispersion  of  flexural  waves  on  lake  ice  agreed  quite  well  with  the  theoret¬ 
ically  derived  dispersion  curves. 

An  expression  for  the  phase  velocity  c  of  flexural  waves  is  easily  derived2  for  the 
case  of  deep  water  (neglecting  air) 


(  1/3)772  72u2  +  (pJpi)gk/Av2y 

1  +  (PjPi  )/2nyy/l  ~  c2/v\~ 


where  y  is  the  ratio  of  ice  thickness  H  to  wavelength  X  f  and  is  related  to  frequency 
/  by  the  equation  y  =  Hf/c  ,  v  is  the  velocity  of  longitudinal  waves  in  the  ice 
sheet,  v*  is  the  speed  of  sound  in  water,  g  is  the  gravitational  constant,  and  Pw,  pi 
are  the  densities  of  water  and  ice  respectively.  It  can  be  seen  from  Eq.(l)  that  for  large 
wavelengths  (y  small),  the  phase  velocity  is  small,  approaching  the  value  ^gX/2  77  of 
gravity  waves  on  deep  water.  As  the  wavelength  decreases  (y  increasing)  the  phase  veloc¬ 
ity  increases.  The  gravity  term  in  Eq.  (1)  can  be  neglected  for  wavelengths  (or  frequen¬ 
cies)  for  which  y  >  0.02.  Group  velocity  V  can  be  computed  as  a  function  of  y  from 
Eq.  (1)  and  the  familiar  expression 

U  =  c  +  y  dc/dy  .  (2) 


Hi* ur ice  Ewing  and  A.  P.  Crary,  "Propagation  of  Elastic  Waves  in  Ice,"  Phyaica,  vol.  5,  no.  7, 
PP  1-10;  1934. 

2  Ibid. 


11 


In  Fig.  2  dimensionless  phase  and  group  velocity  curves  computed  from  Eqs.  (1)  and  (2) 
are  plotted  against  y  for  the  case  vp  =  11,500  ft/sec,  vw  =  4650  ft/sec,  pjp{  ~  1.0904. 
The  dispersion  observed  on  the  seismogram  of  Fig.  la  is  plotted  as  Fig.  3.  Good  agreement 
with  the  theoretical  group  velocity  curves  for  H  =  1.0-1. 1  ft  is  found.  Records  of 
flexural  waves  illustrated  in  Fig.  la  were  obtained  for  all  shots  located  either  in  the 
ice  or  in  the  water.  The  amplitudes  decreased  as  the  depth  of  the  shot  increased,  waves 
of  higher  frequency  decreasing  more  rapidly  than  those  of  lower  frequency,  as  one  would 
expect  from  theory.  Dispersion  was  found  to  be  independent  of  distance  provided  shot 
distances  were  sufficiently  large  for  the  dispersion  pattern  to  develop  (>  500  ft). 

AIR-COUPLED  FLEXURAL  WAVES 


A  pronounced  change  in  record  character  was  observed  for  shots  fired  on  or  above  the 
ice  surface.  Figure  lb  is  the  seismogram  of  a  shot  of  ?  lb  exploded  in  air  17  inches 
above  the  surface  of  the  ice  at  the  same  location  as  the  water  shot  illustrated  in  Fig. 
la.  Both  records  are  aligned  with  respect  to  time  after  the  shot  instant.  Flexural  vi¬ 
brations  from  the  water  shot  are  characterized  by  a  long  dispersive  train  beginning  with 
waves  of  frequency  75  c/s  and  velocity  about  2150  ft/sec  and  ending  gradually  with  waves 
having  a  frequency  of  8  c/s  and  velocity  about  650  ft/sec.  The  air  shot  consists  almost 
entirely  of  a  train  of  constant  frequency  vibrations  the  first  waves  of  which  travel  with 
a  velocity  of  about  2000  ft/sec.  A  gradual  increase  of  amplitude  with  time  is  observed 
(although  some  grouping  tends  to  occur).  The  constant  frequency  waves  build  up  until  the 
arrival  of  a  pulse  (indicated  by  arrows),  after  which  the  amplitudes  diminish  rapidly. 

Figure  4,  a  seismogram,  recorded  at  reduced  gain,  of  an  air  shot  of  2  lb  detonated 
1  ft  above  the  ice  at  a  spread  distance  of  6114-6514  ft,  clearly  shows  the  nature  of  this 
pulse  (indicated  by  arrows).  It  consists  of  a  single  cycle  wave  travelling  with  the  speed 
of  sound  in  air.  There  is  little  activity  in  the  air  shot  after  the  arrival  of  the  air 
wave,  whereas  in  the  water  shot,  flexural  vibrations  continue  to  arrive  after  the  time  cor¬ 
responding  to  that  of  the  air  wave. 

Except  for  minor  changes,  the  general  characteristics  of  the  records  in  Fig.  la  were 
observed  for  all  shots  within  the  ice  or  in  the  water  near  the  ice  bottom.  All  air  shots 
produced  records  similar  to  the  seismograms  of  Figs,  lb  and  4. 

The  constant  frequency  train  of  waves  culminating  with  the  arrival  of  the  air  wave  is 
interpreted  as  an  air-coupled  flexural  wave. 

A  simple  yet  general  explanation  of  the  main  characteristics  of  air-coupled  surface 


12 


waves  can  be  given  with  the  aid  of  the  classical  theory  of  travelling  disturbances.3,4’* * 4 5 
If  the  reaction  of  the  surface  wave  on  the  air  is  neglected,  the  air  wave  can  be  treated 
as  a  pressure  pulse  travelling  over  the  surface  of  a  dispersive  medium  with  velocity  Vg. 
We  can  think  of  the  air  wave  as  imparting  a  succession  of  impulses  at  infinitesimally 
small,  equal  intervals  of  time  as  it  moves  over  the  surface.  Each  impulse  produces  a 
dispersive  train  of  waves.  Superposition  of  the  successive  wave  systems  results  in  rein¬ 
forcement  only  for  those  waves  whose  phase  velocity  is  equal  to  the  velocity  v  of  the 
air  pulse.  The  waves  thus  produced  have  a  discrete  frequency  which  can  be  obtained  from 
the  phase  velocity  curve  of  the  dispersive  medium  from  the  value  y  corresponding  to 
c  =  v  (Fig.  2).  Since  waves  of  a  given  frequency  propagate  with  the  group  velocity  asso¬ 
ciated  with  that  frequency,  the  induced  constant  frequency  surface  wave  will  precede  or 
follow  the  air  wave  according  as  V  >  c  or  V  <  c.  Group  velocity  exceeds  phase  velocity 
for  flexural  waves,  and  accordingly  an  air  shot  can  be  expected  to  produce  a  train  of 
constant  frequency  waves  ending  abruptly  with  the  air  wave.  Air-coupled  gravity  waves  and 
surface  elastic  waves  would  follow  the  arrival  of  the  air  wave. 

Possible  existence  of  air-coupled  Rayleigh  waves  was  discussed  briefly  in  Bateman.6 

If  the  air-ice-water  system  is  treated  as  a  single  acoustic  unit,  one  can  derive  phase 
and  group  velocity  curves  analogous  to  those  presented  in  Fig.  2.  N.  Haskell7  derived 
such  curves  for  the  steady  state.  Identical  curves  were  derived  by  Press  and  Ewing  in  a 
subsequent  investigation  of  the  exact  theory  of  air-coupled  flexural  waves  originating  in 
an  impulsive  point  source  (Part  II,  Fig.  9  (note  that  v  -  t>a)).  A  striking  feature  of 
these  curves  is  the  extreme  steepness  in  the  group  velocity  at  y  .  Since  the  predominant 
waves  observed  at  a  distant  point  are  those  whose  frequency  and  velocity  are  determined 
by  the  group  velocity  curve,8  we  can  expect  a  train  of  waves  of  constant  frequency 
f  a  ~  y av a/H  beginning  at  a  time  r  =  r/2.2vg  and  continuing  until  the  time  t  -  r/v  , 
Calculation  of  the  excitation  function  for  these  waves  (Part  II)  indicates  that  in  an  air 
shot  they  will  predominate  over  waves  whose  travel  time  and  frequency  are  determined  from 


H.  Lamb,  "On  Waves  Due  to  a  Travelling  Disturbance  with  an  Application  to  Waves  in  Superposed  Flu¬ 

ids,"  Phil.  Mag.,  vol.  31,  no.  6,  p  387;  1916. 

4H.  Lamb,  "On  Wave  Patterns  Due  to  a  Travelling  Disturbance,"  Phil.  Mag.,  vol.  31,  no.  6,  p  539' 
1916. 

5  H .  Lamb,  "Hydrodynamics,"  Cambridge  University  Press,  6th  Ed.,  1932  :  pp  413-415. 

6  n  f? 

H.  Bateman,  Rayleigh  Waves,  National  Ac ademy  of  Science  Pr oceedings ,  vol.  24,  pp  315-320;  193 8. 

7 

N.  Haskell,  personal  communication;  March  10,  1950. 

8  fi 

See,  for  example,  T.  H.  Havelock,  The  Propagation  of  Disturbances  in  Dispersive  Media,"  Cambridge 
Tracts  in  Mathematical  Physics,  no.  17,  Cambridge  University  Press;  1914. 
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the  remaining  portions  of  the  group  velocity  curve.  For  an  ice  or  water  shot  the  complete 
theory  indicates  that  the  constant  frequency  waves  are  vanishingly  small  and  waves  follow¬ 
ing  the  part  of  the  group  velocity  curve  to  the  left  of  the  maximum  predominate.  It  should 
be  noted  that  the  dispersion  in  these  waves  for  frequency  /  <  /  is  almost  identical  to 
that  obtained  from  the  group  velocity  curve  of  Fig.  2,  which  was  calculated  under  the  as¬ 
sumption  that  air  had  negligible  effect  on  flexural  waves.  Thus  it  is  only  for  an  air 
shot  that  waves  propagated  according  to  the  steep  portion  of  the  group  velocity  curve 
would  be  observed.  It  is  these  waves  which  we  call  air-coupled  flexural  waves. 

The  frequency  of  air-coupled  flexural  waves  depends  only  on  the  parameter  7a  ,  the 
ice  thickness,  and  the  velocity  of  sound  in  air.  The  value  of  y g  depends  on  the  elastic 
constants  of  the  ice,  in  accordance  with  Eq.  (1).  For  smooth  ice  on  Lake  Superior 
y a  ~  0.091,  and  taking  vg  -  1070  ft/sec,  we  can  write  the  following  simple  equation  for 
frequency  in  terms  of  ice  thickness  H  (measured  in  feet): 

fa  =  91 /H  .  (3) 

Values  of  y g  as  a  function  of  vp  have  been  computed,  assuming  the  same  constant  values 
for  the  other  parameters  of  Eq.  (1)  used  earlier.  The  results  of  these  calculations,  plot¬ 
ted  in  Fig.  5,  can  be  used  to  obtain  y correct  to  within  a  few  percent. 

On  the  smooth  ice  of  Lake  Superior  the  observed  frequency  /  (Figs,  lb  and  3)  averaged 
about  78  c/s.  Ice  thickness  sampled  at  5  locations  along  the  profile  gave  an  average  value 
of  1.1  ft,  and  the  theoretically  expected  frequency  calculated  for  this  value  from  Eq.  (3) 
is  88  c/s,  giving  reasonably  good  agreement,  in  view  of  the  uncertainty  of  the  measured 
value  of  H. 

Records  obtained  on  Lake  Cayuga  gave  a  frequency  of  150  c/s  for  the  air-coupled  flex¬ 
ural  vibrations.  Figure  6  illustrates  a  shot  in  which  a  cap  was  detonated  at  a  depth  of  1 
inch  in  ice  approximately  8  inches  thick.  Water  depth  was  4.3  ft  and  a  spread  of  500-710 
ft  was  used.  Since  -  10,200  ft/sec  on  the  Lake  Cayuga  ice,  yg  ~  0.097  and  we  can 
write  the  following  equation  for  frequency: 

fa  =  m/H  .  (4) 

Since  H  ~  0.66  ft,  /  =  158  c/s.  Agreement  with  the  observed  value  is  good  despite  the 
fact  that  shallow  water  underlay  the  ice  and  that  Eq.  (4)  is  based  on  deep  water  theory. 
The  train  of  dispersive,  long-period,  flexural  waves  which  follows  the  air  wave  on  the  Lake 
Cayuga  air-shot  is  probably  due  to  the  fact  that  the  shot  was  partly  in  the  ice,  thus 
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exciting  both  air-coupled  waves  and  ordinary  flexural  waves. 

Records  of  flexural  waves  from  both  water  and  air  shots  were  obtained  for  propagation 
paths  across  rough  ice  in  which  the  ice  thickness  varied  greatly.  The  flexural  vibrations 
from  a  water  shot  in  rough  ice  exhibited  the  characteristic  dispersion  although  the  waves 
were  more  irregular  and  amplitude  variations  were  more  pronounced.  An  air  shot  across 
rough  ice  produced  the  seismogram,  shown  in  Fig.  7.  In  this  shot,  a  |  lb  charge  was 
detonated  on  the  surface  of  the  ice  at  an  approximate  distance  of  4000  ft  from  the  nearest 
geophone.  A  predominant  frequency  averaging  68  ±  3  c/s  is  observed  for  the  air-coupled 
waves.  Since  the  elastic  constants  of  rough  ice  are  slightly  different  from  those  of 
smooth  ice  the  value  of  y  g  for  the  two  media  will  differ.  For  longitudinal  waves  in 
rough  ice  we  take  the  observed  value  vp  -  10,000  ft/sec.  Using  the  same  values  for 
p  ,  p  ,  pg  ,  ,  v g  as  before,  we  find  y  g  -  0.099,  which  gives  for  the  rough  ice 

fa  =  106  /H  .  (5> 

Thickness  of  the  rough  ice  measured  at  eight  locations  varied  from  1.  1  to  2. 1  ft  and  av¬ 
eraged  1.7  ft.  Computing  the  theoretical  frequency  from  Eq.  (5)  gives  62  c/s,  in  good 
agreement  with  the  observed  value  of  68  c/s. 

Air-coupled  flexural  vibrations  in  rough  ice  show  characteristics  similar  to  those  ob¬ 
served  for  shots  in  smooth  ice.  Wave  amplitudes  build  up  as  the  time  for  the  arrival  of 
the  air  wave  is  approached.  However,  the  air  wave  does  not  arrive  as  a  distinct  pulse  as 
is  the  usual  case  for  smooth  ice  and  some  low  frequency  activity  continues  after  the  air 
wave.  Only  the  microphone  trace,  in  Fig.  7,  registers  the  exact  instant  of  the  passage 
of  the  air  pulse.  Scattering,  air-wave  echoes,  variations  in  group  velocity--all  due  to 
irregularities  in  ice  thickness-- are  believed  to  be  responsible  for  these  differences  on 
the  records  obtained  in  rough  ice.  A  striking  feature  of  the  seismogram  in  Fig.  7  is  the 
recording  of  a  portion  of  the  constant  frequency  wave  train  by  a  microphone  located  ap¬ 
proximately  1  ft  above  the  ice.  Records  have  been  obtained  in  which  microphones  have 
registered  even  larger  portions  of  air-coupled  vibrations.  This  is  indeed  confirmation  of 
coupling  between  ice  and  air. 

As  an  additional  test,  a  j  lb  charge  was  detonated  in  a  tree  on  shore  at  a  distance  of 
20  ft  from  the  edge  of  the  ice  sheet  and  at  a  height  of  20  ft  above  the  ice  surface.  Un¬ 
like  the  usual  air  shot  records  the  seismogram  obtained  showed  no  longitudinal  ice  wave  or 
compressional  water  wave.  However,  the  air-coupled  flexural  waves  were  almost  indistin¬ 
guishable  from  those  obtained  for  shot  points  over  ice. 

In  several  water  shots  for  shot  depths  ranging  from  0  to  50  ft  beneath  the  bottom  of 
the  ice,  a  brief,  faint,  yet  definite  train  of  air-coupled  waves  is  observed  riding  on  the 
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larger  amplitude  flexural  waves.  Since  the  shallow  shots  blow  out,  (i.e. ,  the  gas  bubble 
breaks  out  and  partially  explodes  in  the  air)  a  partial  air  wave  is  set  up,  inducing  a 
coupled  vibration.  For  the  deeper  shots  which  do  not  blow  out  it  is  tentatively  sug¬ 
gested  that  some  sound  energy  enters  the  air  over  the  shot  point  by  diffraction.  Aided  by 
a  favorable  air  velocity  gradient,  this  energy  induces  a  coupled  vibration. 
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PART  II 

THEORY  OF  AIR-COUPLED  FLEXURAL  WAVES 

INTRODUCTION 

In  the  course  of  a  study  of  the  phenomena  associated  with  the  famous  explosion  of  the 
volcano  Krakatoa9,10  in  the  Straits  of  Sunda,  the  authors  noted  a  simultaneous  arrival  of 
an  air  wave  and  tidal  disturbance  at  several  widely  separated  locations.  While  this  curious 
phenomenon  had  been  noticed  previously,11  it  had  been  passed  off  as  a  coincidence;  the  tidal 
waves  were  attributed  to  non-related  earthquakes.  It  is  now  the  belief  of  the  authors  that 
the  simultaneous  arrival  of  an  air  wave  and  tidal  disturbance  resulted  from  coupling  between 
the  atmospheric  pressure  wave  and  the  ocean. 

It  is  surprising  at  first  to  think  of  coupling  between  the  atmosphere  and  the  ocean  es¬ 
pecially  since  the  density  contrast  between  the  two  media  is  so  great.  It  can  be  generally 
shown,  however,  that  coupling  between  air  waves  and  surface  waves  of  all  types  is  ap¬ 
preciable  when  the  phase  velocity  of  the  surface  wave  is  very  close  to  the  speed  of  sound 
in  air. 

The  existence  of  air-coupled  surface  waves  was  firmly  established  in  the  experimental 
studies  reported  in  Part  I  of  this  paper.  A  general  explanation  of  these  waves  was  given 
on  pages  7  and  8  of  that  part. 

In  Part  II  the  exact  theory  for  air-coupled  flexural  waves  in  a  floating  ice  sheet  is 
derived  for  the  case  of  an  impulsive  point  source  situated  either  in  the  air  or  in  the  water 
beneath  the  ice. 

A  complete  discussion  of  the  theory  of  air-coupled  gravity  waves,  with  application  to 
the  explosion  of  Krakatoa,  is  reserved  for  another  paper. 


G.  J.  Symond ,  The  Eruption  of  Krakatoa  and  Subsequent  Phenomena,  TrGbner  &  Co. 

1 0  ft 

C.  L.  Pekeris,  .  The  Propagation  of  a  Pulse  in  the  Atmosphere,"  Proc.  Roy,  Soc 
1  939. 

R.  D.  M.  Verbeek,  Krakatau,  Batavia;  1886. 


London;  1888. 
v  o  1 .  A171,  p  4  34; 
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THE  THEORY 

Consider  the  propagation  of  flexural  waves  in  a  plate  of  infinite  extent  floating  on 
deep  water,  the  thickness  of  the  plate  H  being  small  compared  to  wavelengths  considered. 
Overlying  the  plate  is  an  infinite  atmosphere  having  density  pl ,  and  sound  velocity  v  . 
The  plate  has  volume  density  p2  ,  longitudinal  wave  velocity  vp;  the  water  has  density 
P3  ,  sound  velocity  v3.  A  Cartesian  coordinate  system  is  chosen  with  the  x,  y  axes  in 
the  equilibrium  plane  of  the  plate  and  the  z  axis  vertically  upward.  We  will  use  the 
coordinates  z  and  r  =  ^/x2  +  y2  and  denote  the  corresponding  displacements  by  w  and 
q.  The  subscripts  l,  2,  3  hereafter  refer  to  the  air,  plate,  and  water  respectively. 

SIMPLE  HARMONIC  POINT  SOURCE  IN  AIR 

We  wish  to  determine  the  vertical  displacement  w  of  the  plate  due  to  the  passage  of  a 
system  of  flexural  waves  which  originate  in  a  point  source  of  sound  waves  at  r  =  0  and 
z  -  d.  Particular  interest  is  in  the  solutions  which  predominate  at  large  distances  from 
the  source.  Assuming  simple  harmonic  motion  exp(~  icot)  we  introduce  the  velocity  poten¬ 
tials  <£j  and  <£3  from  which  the  component  velocities  q  and  w  and  the  pressure  p  can 
be  obtained  as  follows: 


p.  -  pfidi.fbt 

q.  =  3<£./3r 

w.  =  ~6c^ifbz 


3  > 


(6) 


It  is  convenient  to  divide  the  air  into  two  regions  by  the  plane  z  -  d  and  to  denote 
values  of  velocity  potential  etc.  for  the  region  0  <  z  <  d  by  primed  symbols. 

It  is  required  that  the  functions  be  solutions  of  the  wave  equation: 


=  32c t./ltt2 


i  =  1 ,  3 


(7) 


V2  =  32/Br2  +  (l/r)B/3r  +  32/Bz2 


where 
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Solutions  of  Eq.  (7)  must  satisfy  the  boundary  conditions  ( for  a  thin  plate) 

3c£' /dz  -  dw2/dt  -  d^/dz  at  z  -  0  (8) 

and  w2  satisfies  the  equation  for  flexural  vibrations  of  the  plate12 

H p2d2w2fdt2  =  -  W3p2v2p  /12)V4u,2  -  p3gw2  ~  p^/dt  +  p^cb'^dt  ,  (9) 

where  v  is  the  velocity  of  longitudinal  waves  in  the  plate,  g  is  the  gravitational  ac¬ 
celeration,  and 


V' 


_LJL 

r  3r 


dr 


n 

r  dr 


r  — 
dr 


Equation  (9)  is  derived  under  the  assumption  that  Poisson’s  constant  for  the  plate  has  the 
value  0. 25. 

We  use  a  method  of  representing  a  simple  harmonic  point  source  originally  given  by 
Lamb. 13  The  procedure  is  first  to  obtain  the  solutions  to  the  problem  where  a  periodic 
pressure  is  applied  to  the  entire  plane  z  -  d  symmetrically  about  the  z  axis  and  then 
to  pass  to  the  case  of  a  point  source  utilizing  the  Fourier-Bessel  integral.  A  point  source 
located  at  r  ~  0,  z  -  d  is  represented  by  ultimately  requiring  continuity  of  pressure  in 
the  plane  z  ~  d  and  continuity  of  vertical  velocity  everywhere  in  the  same  plane  except 
at  the  source  where  the  air  above  and  below  the  source  moves  in  opposite  directions.  Here 
the  discontinuity  in  vertical  velocity  is  proportional  to  a  function  F(r )  which  vanishes 
everywhere  except  at  r  =  0  where  it  becomes  infinite  in  such  a  manner  that  its  integral 
over  the  plane  z  -  d  is  finite. 

Typical  solutions  of  Eq.  (7)  are  of  the  form: 

c£j  =  A  exp(~T) z)JQ(kr)  exp (~icot)  z  >  d  (  10) 

2  2 

Maurice  Ewing  and  A.  P.  Crary,  "Propagation  of  Elastic  Waves  in  Ice,"  Physics,  vol.  5,  part  II, 
no.  6,  pp  1-10;  V934 . 

3H .  Lamb,  "On  the  Propagation  of  Tremors  over  the  Surface  of  an  Elastic  Solid,"  Phil.  Trans.  Roy. 

Soc.  Lond.,  Ser.  A,  vol.  203,  pp  1-42;  1904. 
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d>[  -  [B  exp (tjz)  +  C  exp(-r)  z)]JQ(kr)  exp(-iajt)  0  <  z  <  d  ( H) 


&3  -  E  exp (Cz)J0(kr)  exp(~i cot)  .  z  <  0  (12) 

Flexural  motion  of  the  plate  is  given  by: 

-  DJQ(kr)  exp(~  icL)t)  ,  (  13) 

The  separation  constants  T),  £  obtained  by  substituting  Eqs.  (10),  (11),  and  (12)  in  Eq. 
(7)  are 


rj2  -  k2  ~  co2  /v\  ,  C2  -  k2~co2/v23  ,  (14) 

and  are  defined  as  positive  real  or  positive  imaginary. 

These  solutions  must  satisfy  two  additional  boundary  conditions  at  the  plane  z  ~  d 
where  the  pressure  is  continuous, 

,  (15) 


and  the  vertical  velocity  is  discontinuous, 

Bc£j/3z  -  -  2YJ0(kr)  exp(-  ico  t)  >  (16) 


the  air  above  and  below  the  plane  moving  in  opposite  directions. 

Tbe  expression  for  <£'  is  a  general  form.  and  d3  have  been  chosen  to  decrease 
exponentially  with  distance  from  the  plate  since  we  are  particularly  interested  in  solutions 
for  which  there  is  no  radiation  from  the  plate  into  the  surrounding  media. 

Five  simultaneous  linear  equations  result  when  Eqs.  (10)  through  (13),  the  solutions  of 
Eq.  (7),  are  substituted  in  Eqs.  (8),  (9),  (15),  and  (16),  the  boundary  conditions.  Solv¬ 
ing  for  A,  B,  C,  D,  E  from  these  equations  gives: 


-y 

r]G{k) 


[ g(k )  exp(-r)  d)  +  G(k)  exp(r)d  )] 


(  17) 


4 
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B  -  — [exp(~i7cf)]  (18) 

V 

-  y 

C  E  — — —  [g(k)  exp(~)7d)]  (19) 

V  G(k) 

0  -  — — —  [(p  /p  )r){,  ico  exp(-rjd)]  (20) 

rjG(k) 

E  =  — — —  [(p./p  Wv  exp(-T?cO]  ,  (21) 

vG  (k)  1  2 


where 

G(k)  =  (p3/p2)a2V  +  77C  (Hoi1  -  tfVfeyi2  -  gPj/Oj)  +  (Pj/Pjl&i2?  ,  (22) 

get)  =  (p3/p3)a>2v  +  vl(Ha>2  -H3vy/12  ~  gP3/p2)  ~  (p,/p2)u2t,  .  (23) 

We  now  generalize  the  discontinuity  of  vertical  velocity  in  the  plane  z  ~  d  by  means 
of  the  Fourier-Bessel  integral 


F(r) 


=  I  J0(kr)kdk\  F(k)J  (kk)\dK 
Jo  Jq 


(24) 


F (k)  is  chosen  to  vanish  for  all  but  infinitesimal  values  of  k  in  such  a  manner  that 

F(k)2rrk  dk 

has  a  finite  constant  value.  Thus  taking  Y  -  k  dk  in  the  Eqs.  (10)  through  (13)  and  (17) 
through  (21)  and  integrating  with  respect  to  k  from  0  to  ®  ,  we  obtain  solutions  which 
satisfy  the  appropriate  conditions  at  z  -  0  and  meet  the  additional  requirements  for 
a  point  source,  namely  continuity  of  pressure  in  the  plane  z  -  d  and  continuity  of  vertical 
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particle  velocity  everywhere  in  this  plane  except  at  the  source  r  -  0  where  a  discontinuity 
in  'bcfcfbz  exists  proportional  to 


which  becomes  infinite  in  such  a  way  that  its  integral  over  the  plane  z  -  d  is  finite. 
The  solution  for  a  periodic  point  source  is  therefore 


w„ 


OJ 

-  exp (-icjOt)J  JQ(kr)k  dk  -^2  (p^p^iool  exp(-r)  d)/G(k)  l 


(25) 


Following  the  procedure  of  Lamb, 14  Sezawa, 15  and  Pekeris,16  Eq.  (25)  will  be  evaluated 
by  integration  along  suitable  contours  in  the  complex  plane  Z  -  k  +  im.  let 


E(Z)  -  ~  (px/p2)(Z2  ~  a)2/v~)1/2ioo  exp  [~d(Z2  -  a)2/v2)1/2] 


and  consider  the  two  integrals  whose  sum  can  be  used  to  evaluate  the  integral  in  Eq.  (25): 


and 


H'0(Zr)[E(Z)/G(Z)]Z 


dZ 


(26) 


jH20(Zr)[E(Z)/G(Z)]Z  dZ 


(27) 


*4H.  Lamb,  nOn  the  Propagation  of  Tremors  over  the  Surface  of  an  Elastic  Solid,"  Phil.  Trans.  Roy. 
Soc.  Lond.,  Ser.  A,  vol.  203,  pp  1-42;  1904. 

Sezawa,  "Love  Waves  Generated  from  a  Source  of  a  Certain  Depth,"  Bull.  Earth.  Res.  Inst.  Tokyo, 
vol.  13,  no.  1,  pp  1-17;  1935. 
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C.  L.  Pekeris,  "Theory  of  Propagation  of  Explosive  Sound  in  Shallow  Water 
27,  Propagation  of  Sound  in  The  Ocean;  1948. 


Geol.  So c.  Amer.,  Mem. 


22 


Integrals  (26)  and  (27)  are  taken  around  paths  C  and  C'  in  the  first  and  fourth  quad¬ 
rants  respectively  (Fig.  8).  G(Z)  is  obtained  from  Eq.  (22).  Integral  (26)  has  branch 
points  at  Z  -  op/vx  and  Z  ~  op/ v3  and  a  pole  at  Z  -  kn  which  is  the  root  of  G(Z )  =  0. 
G(kn)  -  0  is  the  characteristic,  or  frequency  equation  for  air-coupled  flexural  waves  on  a 
floating  ice  sheet.  Integral  (27)  has  neither  branch  points  nor  poles.  That  the  singular 
points  lie  in  the  first  quadrant  only  can  be  demonstrated17  by  temporarily  including  small 
frictional  forces  in  the  original  equations  of  motion.  We  denote  by  r\‘  ,  +  £/  the  values 
assumed  by  77  ,  £  on  both  sides  of  the  line  k  ~  op/ v^,  and,  by  ±77"  ,  £"  the  values  on 
the  two  sides  of  k  ~  co/v j  (Fig.  8).  For  small  values  of  m  we  can  write: 

77'  =  1  (cJ/v*  ~  cu2/v2)1/2  >  77"  =  (2moj/vj)1/2exp(  177/4) 


-  ( 2ma>/v3)l/  exp(  177/4) 


i"  =  {co2/v\  ~op2/v\)x/2 


Taking  Integral  (26)  around  the  indicated  contours  and  letting  the  radius  of  the  circular 
arc  approach  infinity,  we  find: 


00  ao 

jH'0(kr)[E(k)/G(k)]k  dk  +fn'0(Z'r)Z"  dn(Pl/p7)a>t," 


exp(~7 fd) 


{(P3/P2)  a>V  +  T,"£'[Wwa-tfVr4/12 -g(p3/p2)]  +  (Pl/p2)a,J£" 


exp(rj"d) 


<p3/p2) o)V-  VC"  (Ho?-  flVz''7i2  ~  gp3/p2)  +  (pJPj)o>2l, 


4  +  f  H'0(Z'r)Z'  dm(Pl/p2)u>  exp(-V'd)  * 


j: 


r 


(P3/P2)wV  +  Tj'CiHo)2 -H3v7Z  V12  -  gp3/p2)  +  (p^oj2^  (p3/p2)co2  T)‘  -  T)'C  (Ho?  -  H*v2Z  4/12  -  gp3/p2)  -  (p\/p2)o)2^' 1 


-/ 


H'0(imr)m  dm(pl/p2)o)  (,n?+  o?/v2)  7  exp[-ld(m2  +o?/v2  )1/J  ] 


i(p3/p2)o>2(m2+o)2/v2)t/2-  (m2+o?/v2)I/1(m2  +  a?/v2)I/2(Ho)2  - H3v2m4 /12  -  gP3/p2)  +  i(p1/p2h?  (n2  +o?/v2)1/ 2 


-  2  TT  i 


tG(k  7/Bfe 


=  0  ,  (28) 


where  kn  is  the  root  of  G(k)~  0  and  Z'  -  o)/v3  +  im  and  Z"  -  op/vx  +  im  respectively 
in  the  second  and  third  integrals  of  Eq.  (28).  Integrating  (27)  around  the  contour  in  the 
fourth  quadrant  gives: 


K.  Sezawa,  "Love  Waves  Generated  from  a  Source  of  a  Certain  Depth,"  Bull.  Earth.  Res.  Inst  Tokyo 
v  o  1  •  13,  no.  1,  pp  1-17;  1935.  ’  ‘  Y  ' 
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< 

/ 


„  Effe) 

H \(kr) -  k  dk 

0  G(k ) 


o 

I 


H^(mr  )m  dm(Px  /p2  )  co  ( m2  +  <v2  /v2  ) 1/2  exp[~  idfm2  +  ai2/t>2  ) 1/2] 


i(p3/p2  )  6J  2(m2  +  w2/v2)l/2  -  (m2  +  a>2/t’2)1/2(m2  +  w2/r2)I/2(/y<x>2- W3v2m4/12  -  gf^/p2)  +  i  )  w2(  m2  +  o)  2/r2  )1/: 


=  0  .  (29) 


The  integrals  along  the  infinite  arcs  in  the  first  and  fourth  quadrants  vanish  because  of 
the  well  known  properties  of  the  Hankel  functions.  From  the  relation  H*(im)  =  -  H2(~im) 

it  can  be  seen  that  the  last  integrals  of  Eqs.  (28)  and  (29)  are  equal  and  opposite  in  sign. 
Since  2 JQ(kr)  ~  Hl0(kr)  +  H^fkr),  the  sum  of  the  first  integrals  in  Eqs.  (28)  and  (29) 
is  equivalent  to  Eq.  (25).  Using  the  expansions  Hq(x)  exp  [  i  fx  -  77/4)] /(t7x/2)1/2  and 
H2(x)  exp[-ifx~  7r/4)]/(77x/2)1/2  ,  we  add  Eqs.  (28)  and  (29)  and  obtain  for  large 

values  of  r 

u>2  -  wQ  +  w'  +  w"  ,  (30) 


where  wQ  represents  the  residue: 


*0 


277  i 


k  Efk  )  exp [i(k  r~  cot  -  77/4)] 

n _ n _ n _ 

(  nkn  r/2)1/2  BG  (kn)  /  d  k 


(31) 


and  w'  and  w"  are  the  two  branch  line  integrals: 


w 


-2  exp[t(roj/i;^  -  cot  -77/4)]  C  exp(~  mr) 


(77/2  )1/: 


rp/2  z'  dm(px/p2)a)  exp(~  T)' d)  * 


(p3/p2)^2v'V 


UpJp.W-n'y  -  -//VZ'4/12  -  tPjpp  +  (pt/pt)  ? 


(32) 


II 


w 


2  exp[i(rco/v1  ~  cot 

(77/2)1/2 


exp(“  mr) 
(Z"  r)1/2 


Z"dm(px/p2)coV 


x 


cc&(n"d)l(pt/p,)0,*v"  +  -  «VZ"Vl2  -  epa/p2)]  + 

-  Up,/pt)v*ri‘  +  7 fV'fHco2  -  HV*"V12-  ept/p2P  + 


sinh(Tj" d)(px/p2  )co2^" 
Up1/p2)u21"]2 
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where  Z'  ~  a)/u3  +  i m  in  Eq.  (32)  and  Z‘  ~  co/ v 2  +  im  in  Eq.  (33).  These  definite 
integrals  can  be  expanded  in  asymptotic  forms  by  means  of  the  following  formula  which  is 
readily  derived  from  the  integral  representation  of  the  gamma  function: 


(m)  exp  (~mr)  dm 


r<V2),,,,  n>  +  rx  5/2i  . 

r3/2  r  5/2  l! 


For  large  r  we  have: 


(34) 


2i(px/o3)  exp  [i(rco/v3  -  co  t )]  exp  [~id(  cjo2 /v2  ~  a?/v2)l/2] 

r 2  co(v2/v2  ~  1)I/2 


(35) 


2  exp[ ifrcu /v}  -  co  t)  ] 


r2v. 


(  *>7/v]  -  cJ/v\)x/ 2 


+  (p7/ox)  [H-H\2cJ/\2v\  -  g(P3/p2)/co2]  +  d 


(36) 


The  solutions  given  by  Eqs.  (30)  through  (36)  and  (35)  through  (36)  have  been  expressed 
as  a  residue  and  two  integrals  along  branch  lines.  The  branch  line  integrals  represent 
waves  which  travel  with  velocities  t>3  and  Vj  and  diminish  rapidly  for  large  ranges  as 
r*2,  The  residue  represents  a  train  of  dispersive  waves,  whose  frequency  equation  is  given 
by: 

G(kn)  =  0  (37) 


Since  these  waves  diminish  in  amplitude  as  r' 1/2  they  represent  the  predominant  disturbance 
at  large  ranges.  Using  the  expression  for  phase  velocity  c  -  co/k  we  can  rewrite  Eq.  (31) 
in  a  form  more  convenient  for  computation: 

W(kn) 

w0 - ri/2v  wi;rexP[“  dkn ^  “  c2/v]?/2]  exp[i(knr  ~  cot  -tt/A)]  ,  (38) 

p‘ 

where 


W(k) 


2(27T)1/Vp1/^2)  (Un,/2(c/viiH\  -  cV^)1/; 


(39) 
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1 

I 


(p2!p2  )c2/ v 


(i  - 


P  + 


(1  -  c2/v\)^2 

(1  -  r2/,,2  )l/2 


(1  ~  c2/ v? 


/v\yn  (i  -  c2/v21)1''2 


kHc2/v 2 

'  p 


~  (kH)3/ 12 


“  (P3/P2)(g/k)/v] 


-  (1  -  c2/v?)‘/2(l  -  c2/v2)'S2(kH)3/3  +  /Pl/P.2.)c-/-L 

1  3  (1  -  c2/vl),/2 


(40) 


Equation  (37)  likewise  may  be  rewritten  as: 


(P2/p,)(c2/v2p)(\  -  c2/v\)'/2  +  (1  -c7up1/2(l  -  cI/i.2V/2[it«c1A2  -  (kH)’/12  -  (p,/p,)(g/k)/ v2 


+  (pjp2)(c2  !x,\)(\  -  <r2/t,2)"2  =  0  .  (41) 


Equation  (41)  is  a  cubic  in  kH  with  the  elastic  constants  of  the  system  as  param¬ 
eters.  It  defines  an  implicit  relationship  between  the  frequency  /  =  ck/2rr  and  the  phase 
velocity  c  of  waves  given  by  Eq.  (38).  For  c  <  <  v3,  real  values  of  kH  exist. 
These  are  plotted  in  Fig.  9  as  a  function  of  the  dimensionless  parameter  y  -  kH/2n  for 
the  case  vt  =  1070  ft/sec,  vp  -  11,500  ft/sec,  v3  -  4650  ft/sec,  P\/P2  ~  0.00141, 
P3/P2  ~  1.0904.  The  factor  Y/(k )  shows  the  dependence  of  amplitude  upon  frequency.  The 
exponential  term  in  Eq.  (38)  shows  the  influence  of  the  height  of  the  source  on  anplitudes. 
It  is  seen  that  wave  amplitudes  decrease  exponentially  as  -  dk  (1  -  c2/v2)1/2 


GENERALIZATION  FOR  AN  ARBITRARY  PULSE 

In  a  dispersive  system,  such  as  the  one  under  discussion,  energy  associated  with  each 
wavelength  or  frequency  is  known  to  propagate  with  the  group  velocity  given  by  the  familiar 
expression: 

U  ~  c  +  k  dc/dk  .  (42) 

Carrying  out  the  indicated  differentiation,  using  the  functional  relationship  between  c 
and  kH  given  by  Eq.  (41),  enables  one  to  compute  values  of  If  which  are  plotted  in 
dimensionless  form  in  Fig.  9. 

Having  obtained  the  steady  state  solution  for  a  simple  harmonic  point  source  in  the 
air,  we  proceed  to  the  case  of  an  arbitrary  initial  disturbance.  If  the  time  variation  of 
the  initial  disturbance  at  the  source  is  f(t )  represented  by  its  Fourier  Transform 
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f(t)  =  (1/27T )  J  exp(-  icot)g(Ci))  d 
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and  the  displacement  w.  becomes: 


OD 

wQ  =  r~  1/2W*  1/2(2rrvp)~  lf  g(oo)  W  (kj  exp[-  dkj  1  -  c2/v2)l/ 2]  exp[i(knr~ut  -  n  /  \)]  doj  .  (45) 

-00 

In  Eq.  (45),  kr  is  a  function  of  oo,  through  Eq.  (41).  We  follow  Pekeris  in  represent¬ 
ing  the  initial  disturbance  created  in  an  explosion  as  follows: 


f(t )  -  expf~  crt)  t  >  0 

=  0  .  t  <  0 


(46) 


With  this  definition  we  can  write 


f(t) 


00 

(1/2  tt)  I  (cr  -  ico)'1  exp(-  icot)  dco 
~  -00 


g(co)  -  (cr  -  ico)'1,  (47) 


(cr  -  io. >  )‘  lW(k  )  exp[~  dk  ( 1  -  c2/v2)1/2]  exp[iffe  r  -  co  t  -  tt/4)]  dcd  •  (48) 

n  n  l  n 

-CD 

To  evaluate  this  integral  we  use  Kelvin’s  approximate  method  of  stationary  phase.  For  a 
discussion  of  this  method,  and  the  validity  of  the  approximations  used  the  reader  is  re¬ 
ferred  to  papers  by  Havelock, 18  Pekeris,19  and  Eckart.20  Integrating  Eq.  (48)  by  the  method 
of  stationary  phase  gives: 


W(kn)  exp [~  cfen(  1  ~  c2/v2)i/2]  exp[iffenr  -  eont  -tt/\  ±7t/4)] 

rHfvJv x)(a-  ia>0 )[( vjc)(vj\})2 d(U/v1  )/dy] 1/2 


(49) 
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where  the  upper  or  lower  sign  in  the  exponential  term  is  to  be  taken  according  as 
d(li/vl)/dy  is  negative  or  positive,  and  kQ,  co0,  U/v1  d(U/vl)/dy  are  evaluated  from 
the  phase  and  group  curve  for  values  of  r  and  t  which  satisfy  r/t  -  U .  The  approxima¬ 
tion  used  to  derive  Eq.  (49)  is  valid  provided  r  is  large  and  t  is  sufficiently  removed 
from  stationary  values  of  group  velocity.  To  complete  the  solution  we  add  to  Eq.  (49)  its 
complex  conjugate  since  there  are  two  stationary  points  ±cu0  (a  solution  to  Eq.  (48)  except 
for  a  reversal  of  sign  of  the  phase  of  the  exponential  term  would  have  been  obtained  had  we 
assumed  initially  a  time  factor  exp(ioit)  instead  of  exp (~icot)).  Thus  we  obtain: 


2W(k  )  exv[~  dk  (l  ~  c2 /v2  )1/ 2] cos (k  r  ~  cot  +  tan'1  co/a) 
w  —  - u - v - — — - y — - — . 

rH  (vp/vx  )(a2  +  a?)1/2  [(v1/c)(v1/U)2d(U/v1  )/dy] 1/2 

for  d(V/v1  )/dy  <  0 

2lf(kn)  exp[~  dkn(  1  ~  c2/t>2  )1/2]  sinffe0r  -  cot  +  tan  ~1(o/a) 
0  rH(vJvl)(a2  +  co2)1/2  [(v1/c)(v1/U  )2d(U/v1 )  jd  y] 1/2 


(50) 


(51) 


for  d(U/vl)dy  >  0  . 

A  glance  at  the  phase  and  group  velocity  curve  in  Fig.  9  shows  that  the  approximation 
used  in  evaluating  Eq.  (48)  fails  at  the  times  tm  -  r/Umgx  and  t  -  r/vx  corresponding 
respectively  to  the  arrival  of  waves  travelling  with  a  maximum  value  of  group  velocity  and 
waves  for  which  y  >  y  g,  where  the  phase  and  group  velocity  are  constant. 

For  the  maximum  value  of  group  velocity  the  method  of  stationary  phase  can  still  be 
used  to  evaluate  Eq.  (48)  and  gives  (see  Havelock21): 

(4/3)sinf^/3)r(l/3)H/ffe0)exp[-(ife0(l  -  c2/v2  )1/2]cosf/eQr  -  coQt  ~  7t/4  +  tan'1  co/a)  (52; 

11,0  rs/6(l/6)1/3ft>p/Vj  )H1/2(a2  +  co2  )1/2[2ttH  2(vx/c)2  (vl/V)2d2  (U/v^  )/dy2]1/3 

In  the  region  y g  <  y  <  00  of  constant  phase  and  group  velocity  the  method  of  stationary 
phase  fails  to  evaluate  the  integral  in  Eq.  (48).  Waves  propagated  at  a  constant  value  of 
phase  and  group  velocity  are  nondispersive,  and  the  signal  arriving  at  a  time  t  -  r/v j 
would  have  the  analogous  characteristics  of  a  pulse  transmitted  through  a  high  pass  filter 
with  a  low  frequency  cutoff  at  f g  -  ygvl/H  .  However,  the  excitation  function  W(k)  is 
vanishingly  small  for  all  /  >  f g,  and  has  appreciable  value  only  at  /  -  /a*  Consequent¬ 
ly,  the  resultant  signal  at  t  ~  r/vx  would  have  the  characteristics  of  a  pulse  transmit¬ 
ted  through  a  sharply  tuned  band  pass  filter  with  peak  frequency  at  f  -  f g. 
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DISCUSSION  OF  THEORY 


In  the  previous  section  the  steady  state  solution,  Eq.  (25),  was  expressed  as  the  sum 
of  the  residue  of  the  integrand  and  two  integrals  along  branch  lines  corresponding  to  branch 
points  k  -  Cl>/v1  and  k  -  oj/v3.  The  branch  line  integrals  represent  waves  travelling  with 
velocities  v3  and  v1  respectively  and  diminishing  in  amplitude  with  distance  as  r~ 2 . 
The  residue  contributes  dispersive  waves  which  predominate  at  large  distance  since  they 
diminish  only  as  r'1/2.  Equation  (36)  gives  the  functional  relationship  between  phase 
velocity  and  frequency  which  characterizes  the  dispersion.  Phase  velocity  is  plotted  as  a 
function  of  the  dimensionless  parameter  y  ~  kH/2n  -  Hf/c  in  Fig.  9.  Wave  amplitudes  for 
the  steady  state  depend  on  frequency  through  the  factor  W(k )  given  in  Eqs.  (39)  through 
(40).  A  graph  of  W(k)  as  a  function  of  y  is  presented  in  Fig.  10,  for  a  source  in  the 
air.  Interchanging  the  subscripts  l  and  3  in  Eqs.  (39)  through  (40)  enables  one  to  compute 
the  corresponding  amplitudes  for  a  source  in  the  water  (Fig.  11). 

Study  of  these  curves  for  the  steady  state  reveals  that  a  peak  amplitude  exists  for  a 
point  source  in  the  air  at  a  frequency  f a  corresponding  to  the  phase  velocity  c  ~  v }  . 
For  a  point  source  in  the  water  largest  amplitudes  are  associated  with  low  frequency  waves. 
As  the  frequency  (and  phase  velocity)  increases,  wave  amplitudes  decrease  and  abruptly  drop 
to  zero  as  the  frequency  /  is  approached. 

Generalization  of  the  steady  state  solutions  for  an  impulsive  point  source  leads  to  the 
group  velocity  curve  presented  in  Fig.  9.  At  a  distant  point  waves  of  a  given  frequency 
will  arrive  at  a  time  corresponding  to  propagation  at  the  group  velocity  associated  with 
that  frequency.  In  addition  wave  amplitudes  are  proportional  to  the  inverse  square  root 
of  the  slope  of  the  group  velocity  curve.  Referring  to  Eqs.  (50)  through  (51),  one  can  ob¬ 
tain  an  amplitude  factor  for  an  impulsive  point  source 


P(K) 


2W(k„)  exp[~  dkjl  -  cVv*)172] 
[(v1/c)(vl/V)2d(U/v1  )/dy]1/2 


(53) 


from  which  wave  amplitudes  at  large  range  can  be  obtained.  Graphs  of  P(kQ)  have  been 
computed  from  Eq.  (48)  and  are  presented  in  Figs.  12  and  13  for  an  air  shot  and  water  shot 
respectively,  located  at  a  distance  d  ~  5H  from  the  ice  sheet.  With  the  aid  of  Figs.  9, 
12,  and  13,  one  can  describe  in  a  fairly  complete  manner  the  sequence  of  waves  arriving  at 
a  distant  point. 

For  an  air  shot  the  first  waves  to  arrive  appear  at  the  time  t  —  r/2.2v1  corresponding 
to  propagation  at  the  maximum  value  of  group  velocity.  These  waves  appear  with  large  ampli¬ 
tudes  and  with  a  frequency  close  to  /  .  As  time  progresses  two  wave  trains  arrive  simul¬ 
taneously  corresponding  to  the  two  branches  of  the  group  velocity  curve  on  either  side  of 
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the  maximum.  Waves  propagated  according  to  the  left  branch  of  the  group  velocity  curve  are 
dispersive  and  undergo  a  rapid  reduction  in  amplitude  as  the  frequency  decreases  from  the 
value  /  .  Waves  propagated  according  to  the  steep  branch  of  the  group  velocity  curve  to 
the  right  of  the  maximum  appear  as  a  constant  frequency  train  continuing  to  the  time 
t  “  r/Vj.  The  constant  frequency  waves  would  appear  with  relatively  large  amplitudes  since 
their  frequency  fg  is  close  to  the  peak  amplitude  frequency.  Thus  the  predominant  dis¬ 
turbance  appearing  in  an  air  shot  would  consist  of  a  train  of  constant  frequency  waves  be¬ 
ginning  at  the  time  t  -  r/2.2vx  and  ending  abruptly  at  the  time  t  =  r/Vj.  At  this  time 
an  air  wave  would  arrive  corresponding  to  the  flat  portion  of  the  group  velocity  curve.  For 
reasons  given  earlier,  it  would  have  the  characteristics  of  a  pulse  transmitted  through  a 
sharply  tuned  band  pass  filter  with  peak  frequency  at  /  . 

For  the  case  of  a  water  shot  the  sequence  of  arrivals  would  be  the  same  since  the  group 
velocity  curve  in  Fig.  9  is  still  applicable.  However,  the  amplitudes  of  the  waves  differ 
greatly  (Fig.  13)  in  that  the  dispersive  waves  corresponding  to  the  left  branch  of  the 
group  velocity  curve  are  predominant.  A  water  shot  would  therefore  be  characterized  by  a 
train  of  dispersive  waves  beginning  at  the  time  t  =  r/2.2v1  with  a  frequency  slightly 
less  than  /  .  As  time  increases,  the  frequency  decreases  and  the  amplitude  increases.  Al¬ 
though  the  steep  portion  of  the  group  velocity  curve  still  contributes  constant  frequency 
waves,  these  are  now  negligibly  small.  It  is  interesting  to  note  that  the  portion  of  the 
group  velocity  curve  to  the  left  of  the  maximum  is  identical  to  that  which  would  have  been 
obtained  had  we  neglected  the  air  and  it  is  this  portion  of  the  group  velocity  curve  which 
is  responsible  for  the  predominant  waves  appearing  in  a  water  shot.  Thus  for  frequencies 
less  than  /  the  dispersion  and  amplitudes  of  flexural  waves  from  a  water  shot  are  unaf¬ 
fected  by  the  air.  It  is  only  for  an  air  shot  that  a  constant  frequency  train  of  "air- 
coupled"  flexural  waves  appears. 
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APPENDIX 

APPLICATIONS  OF  AIR-COUPLED  FLEXURAL  WAVES 


Air- coupled  flexural  waves  can  provide  a  simple  means  for  determining  the  thickness  of 
a  floating  ice  sheet.  If  the  elastic  constants  of  the  different  types  of  floating  ice 
(i.e.,  smooth  ice,  rough  ice,  lake  ice,  salt  water  ice,  etc.)  in  a  given  area  are  known,  one 
can  obtain  the  value  y from  a  curve  such  as  that  given  in  Fig.  5.  This  leads  to  a  simple 
numerical  relationship  between  frequency  and  average  ice  thickness  between  shot  point  and 
detector.  If  the  elastic  constants  are  unknown,  the  observed  value  of  f  can  be  used  to 
estimate  the  flexural  rigidity  of  the  ice  sheet.  Ordinary  caps  can  provide  sufficient 
energy  to  obtain  average  ice  thickness  over  profiles  of  a  few  hundred  feet  to  several  thou¬ 
sand  feet.  Larger  charges  may  be  needed  for  larger  distances. 

Air- coupled  flexural  waves  are  particularly  adaptable  for  use  with  radio  monitored 
geophones  parachuted  from  a  plane. 

Additional  tests  are  needed  to  gauge  the  accuracy  of  measurements  of  ice  thickness  using 
air- coupled  flexural  vibrations.  Present  data  suggest  that  thickness  can  be  determined  to 
within  10  percent,  possibly  better. 
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SECONDS  25  corr.  ♦  .061  sec. 
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Fig.  2.  Flexural  waves  from  a  water  shot  (a")  and  air-coupled  flexural  waves  from  an  air  shot 
(b)  in  smooth  ice  1.1  ft  thick. 
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Fig.  3.  Observed  and  theoretical  dispersion  of  flexural  waves 
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Fig.  4.  Air-coupled  flexural  waves  in  smooth  ice  1.1  ft  thick  recorded  at  reduced  gain. 
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Fig.  5.  Frequency  parameter  y a  plotted  as  a  function  of  longitudinal  wave  velocity. 
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Fig.  6.  Air-coupled  flexural  waves  in  smooth  ice  .66  ft  thick. 


Fig.  7.  Air-coupled  flexural  waves  in  rough  ice  averaging  1.7  ft  in  thickness. 
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Fig.  8.  Integration  paths  in  the  complex  k  plane. 
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Fig.  9.  Dimensionless  phase  velocity  (c/vl)  and  group  velocity  ( ll/v 
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Fig.  10.  Steady  state  amplitude  function  for  an  air  source. 
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Fig.  11.  Steady  state  amplitude  function  for  a  water  source. 
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Fig.  12.  Amplitude  function  for  an  air  shot  at 
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Fig.  13.  Amplitude  function  for  a  water  shot  at 
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